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O , It is commonly accepted that quasar redshifts have a cosmological character and 

, that most of the quasars are at Gigaparsec distances. However, there are some cases 

\ where several quasars with completely different redshifts and a nearby active galaxy are 

\ aligned in a certain way or occupy a very small patch on the sky, which is claimed by 

^ I some authors to be unlikely to happen by chance. Is there a small subset of quasars 

' with non-cosmological redshifts? For quasars apparently associated with galaxies, we 

■ consider two scenarios for the origin of their redshift: 1. a change in the scale factor 
lO ■ 

Q>^ . of the whole Universe (standard, cosmological scenario), 2. a velocity-induced Doppler 

, shift of a nearby object's spectrum (local, ejection scenario). We argue for a simple 

i astrometric test which can distinguish between these two sources of quasar redshifts by 

"fH \ constraining their proper motions. 

Qh| We give the predictions for the maximum possible proper motions of a quasar for 

O ' the cosmological and local scenarios of the origin of their redshifts. We apply these 

' theoretical results to the Bukhmastova (2001) catalog, which contains more than 8000 

■ close QSO-galaxy associations. In the standard interpretation of quasar redshifts, their 
. typical proper motions are a fraction of micro arc-second, and beyond the reach of 

^ , planned astrometric missions like GAIA and SIM. On the other hand, the quasars 

^ \ ejected from local AGNs at velocities close to the speed of light would have proper 

■ " " ' motions 5-6 orders of magnitude larger, which would easily be measurable with future 

astrometric missions, or even in some cases with HST, VLT and Keck telescope. The 
distributions of proper motions for the cosmological and local scenarios are very well 
separated. Moreover, the division corresponds nicely to the expected accuracy from 
GAIA and SIM. 

Subject Headings: astrometry — cosmology: miscellaneous — galaxies: active — 
quasars: general — relativity 



1. Introduction 



Quasi-stellar objects (QSO, quasars) were identified as a new class by Schmidt (1963). The 
first object with an understood spectrum was 3C 273, with a redshift of 0.158. Many more followed 
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and very soon (Burbidge k, Burbidge 1967) it was obvious that there are no quasars with observed 
blueshifts. A wavelength shift in quasars can be attributed to: 

1. a cosmological redshift, i.e. due to expansion of the whole Universe (hereafter referred to as 
the cosmological scenario), 

2. a gravitational redshift due to high concentration of material in the emission region, 

3. a pure Doppler shift associated with the rapid relative motions of the system quasar-observer 
(hereafter referred to as the local scenario), 

4. some, as of yet undiscovered, phenomenon. 

Already in XIV century William Ockham realized that it is not necessary to create new entities 
if the explanatory potential of the known entities is sufficient — therefore, we will exclude the 
last possibility from the current discussion. At the first glance, the absence of blueshifted quasars 
strongly suggests that quasar redshifts have either a cosmological or gravitational origin. The 
conditions needed to produce gravitational redshifts are so extreme that we do not discuss this 
possibility any further. Historically, a cosmological character of redshifts was criticized based on 
the argument that maintaining the observed quasar fluxes requires huge energy supplies. This issue 
is no longer perceived as a problem, since the accretion of material on a central, super-massive black 
hole provides a satisfactory explanation. The Doppler (ejection) origin for quasars redshifts would 
be a viable alternative only if one could devise a way to avoid the blueshift problem. In the 
presence of only the longitudinal Doppler effect, about half of the randomly-ejected objects should 
be blueshifted. The situation, however, changes dramatically when the transverse Doppler effect 
becomes important. In this case, many of the objects approaching the observer are redshifted rather 
than blueshifted. Still, in a significant population of quasars, one should see a few objects that are 
blueshifted. 

Current observational techniques have enabled well defined selection processes for QSO's. Sev- 
eral large area sky surveys like 2 Degree Field and Sloan Digital Sky Survey have opened a new 
window in the identification of QSOs and allowed the acquisition of a wealth of information (Croom 
et al. 2001; Schneider et al. 2002). All of the thousands of quasars that we know of are redshifted. It 
is unrealistic to expect that the ejection process would always meet the conditions to produce only 
redshifts. Therefore, a universal explanation of the redshifts of all quasars by the local hypothesis 
would require a giant conspiracy^. But the ejection mechanism may be responsible for the redshifts 
of a subset of QSOs. There arc some cases where several quasars with completely different redshifts 



^For individual quasars, the additional evidence for their cosmological redshift may come from the identification 

of their host galaxies in deep imaging or the detection of Lyman-a forest in their spectra. However, from the ground, 
the Lyman-a detections are possible only for the quasars with a> 1.6 (e.g., Kim et al. 2003), which is due to the 
cutoff associated with atmospheric absorption. 
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and a nearby active galaxy occupy a very small patch on the sky, which is claimed to be unlikely 
to happen by chance (e.g., Arp 1999). Also, it is suggested that X-ray compact sources tend to 
group around nearby active Seyfert galaxies (Burbidge, Burbidge, &: Arp 2003). It is even claimed 
by some that they preferentially lie along the minor axis of a Seyfert, and therefore were probably 
ejected from it (Chu et al. 1998, Arp 1999). Here we do not engage in another discussion about the 
significance of such cases. However, we also do not follow the route of the overwhelming majority of 
the astronomical community which ignores this evidence. Instead wc entertain the idea of a simple 
astrometric test which can produce a clear, conclusive answer. Wc show that measuring the proper 
motion (p) of QSOs or at least placing an upper limit on /j, will allow one to distinguish between 
Doppler kinematic redshifts and the cosmological ones. 

Before presenting more detailed computations, it is instructive to see why the measurements of 
quasar proper motions promise to distinguish between the local and cosmological origin of quasar 
redshifts. In the standard cold dark matter (CDM) cosmology the peculiar velocities of high- 
redshift galaxies, are expected to be only a tiny perturbation on top of the overwhelming Hubble 
flow. The same is true for quasars because they form in the centers of galaxies. Therefore, the 
peculiar velocities of quasars in the cosmological scenario are likely to be ug^cosmo ~ lO^'^c, where c 
is the speed of light. On the other hand, if the local ejections are to produce redshifts, the ejection 
velocities must be close to the speed of light, i.e. f (5,iocai ~ c. In addition, the local quasars would be 
at a distance of -DQ,iocai ~ 10 Mpc, and the cosmological ones at Dq 

cosmo 1 Gpc. Consequently: 

/^QjCosmo ^ -Dq^cositio -J — 1 |-'^Q,cosmo 10~^ (1) 
A*Q,local -CQ,local ^-'Q, local 

The idea of the astrometric test we present exploits this large difference between the expected 
proper motions. In §2 we briefly review the theory of relativistic ejection and derive some new 
results. Future astrometric measurements that will allow the implementation of the proper motion 
test are discussed in §3. In §4 we describe the catalog of QSO-galaxy associations that we use, and 
in §5 derive the expected proper motions of quasars in this catalog. We summarize the results in 
§6. 



2. Theoretical considerations 

We start with a brief review of the kinematics of relativistic ejection (Bchr ct al. 1976). We 
note here that the mathematical formalism of relativistic ejection can be used to describe both the 
local and the cosmological case. In the local scenario, a quasar is ejected with a relativistic velocity 
from an active galaxy. In the cosmological scenario, a quasar may be viewed as being "ejected" 
with its small peculiar velocity from the frame of reference moving with the Hubble flow. 

We introduce a cosmological frame of reference with the observer at rest at its origin. The 
line connecting the observer and some point G (which may or may not represent a galaxy) , defines 
the line of sight. A quasar Q is ejected from point G with a velocity represented by /3 = v/c as 
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measured with respect to the local cosmological standard of rest at G. The angle of ejection a in 
the G rest frame is the angle between the quasar velocity and the direction along the line of sight 
away from the observer. With these definitions one obtains for the proper motion /j, oi Q with 
respect to G, 

c /3sina , , 

^~ 'D^1 + Pcosa' ^ ^ 

where Dm is the transverse comoving distance or proper motion distance. At this point we choose 
standard A cosmology with ^Im = 0.3 and Q\ = 0.7 (e.g., Bennett et al. 2003), and set the curvature 
contribution to the total mass density 0,k = 0. Then the transverse comoving distance is defined 
as 

where the distance Djj = c/Hq is the Hubble radius^, and we take the Hubble constant to be 
Hq = hx 100 km/s/Mpc. Under our assumptions, the function E{z) (e.g., Hogg 2000) is given by: 

E{z) = ^nM{i + z)^ + nA. (4) 

From (2) we see that the dependence of the proper motion on the local geometry of ejection 
process is described by the function F{a,f3). 

(3 sin a , , 

1 + p cos a 

The proper motion in the cosmological scenario can be derived from Eq. (2), if we associate ejection 
velocity with the peculiar velocity of the quasar. In this case, we simply consider "an ejection" from 
the reference frame moving with the Hubble flow. For a given P = Pq a maximal proper motion is 
achieved for the angle ao such that cos ao = Po and is equal to 

Mmax ~ n o 1 /2 

Dm (1-/3o2)'/' 

We assume that the peculiar velocity does not exceed (3o = (1000km/s)/c. This is very likely to be 
a conservative assumption because the pairwise velocity dispersions at 2; ~ are of the order of a 
few hundred km/s (Jing &; Borner 2001; Landy 2002) and they do not rise with z for the currently 
favored A cosmology (Barrow & Saich 1993). Hence, /icosmo is very likely smaller than: 

/^cosmOjinax — 0.70367 X 10"^ X — ^ arcsec/yr. (7) 

Dm 



^Note that in equation (3) the distance Dm is expressed in terms of the redshift z of point G, which we will 
approximate by the measured redshift of a quasar in the cosmological scenario or by the redshift of a galaxy in the 
local scenario. This approximation leads to negligible errors in the cosmological scenario, but may be important for 
local galaxies, since they have peculiar radial velocities that are often a substantial fraction of the radial velocity of 
point G with respect to the observer. In some instances, when the appropriate observational data exist, this problem 
may be remedied by using measured distances to galaxies instead of the ones derived from their redshifts. 
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The fraction Dm/Dh as a function of redshift for A cosmology with Qm = 0.3 and Ua = 0.7 is 
displayed in Figure 1. 

In the local scenario the ejection is caused by a host galaxy. For an ejected QSO (subscript 
Q) and the galaxy (subscript G), Behr et al. (1976) define: 

1 + ZQ l + /3cosa 

K=- - = —, (8) 

1 + ZG (1-/32)^ 

where zq and zq are the redshifts of the QSO and galaxy, respectively. Since a typical quasar has 
ZQ ~ 2, and a typical galaxy with which it is apparently associated has zq <C 1, the coefhcient k 
peaks around 3. 

Our goal now is to derive the maximum possible proper motion given the angle of ejection a 
and coefficient k. Eq. (8) yields 



_ cos a + H V /V- - sill- a 
+ cos^ a 

Substitution of eq. (9) into (5) yields 

, sinacosa + Ksinav — sin^ K , , 

F{a,K) = -. (10) 

+ 2 cos^ a + cos any — sin^ k 

We have plotted F{a, k) as a function of a for k = 2, 3 and 4 (Fig. 2). For each k, F{a, k) peaks 
at a characteristic angle ttmax given by 

- 1 arccos[-^^^^] if . > ^2. ^''^ 

The algebraic manipulation of Eqs. (10) and (11) leads to the conclusion that the maximum 
observable proper motion in the local scenario is given by: 

C K / -t n\ 

Wocal.max 7:; 77 > V / 

which is a very interesting new result. 



3. Astrometric Missions 

The aim of this analysis is to design a test that will use the proper motions to interpret 
quasar redshifts. Therefore, one should consider the most powerful facilities that can make such a 
measurement. We limit our investigation to the upcoming astrometric missions GAIA^ and SIM^ 



^http : // astro . estec . esa . nl/GAIA/ 

*http : //planetquest . jpl .nasa . gov/ SIM/sim_ index . html 
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that are supposed to be launched around 2010. Both will reach a parallax (positional) measurement 
accuracy of about a few micro arc-seconds (;Uas) at visual magnitude V = 15. SIM will be able to 
point at tens of thousands of stars and GAIA will scan millions of stars over the entire sky. Here 
we describe the GAIA accuracy in some detail. We assume that SIM can reach similar accuracy 
with the proper integration time; moreover we expect that integration time that would allow SIM 
to surpass GAIA substantially will in practice be very hard to obtain. The parallax accuracy as a 
function of magnitude m for a mission with a GAIA-like observational setup is given by Gould &i 
Salim (2002): 

CTpa. = ^C7o2 + fTig^ X 100-4(--15) [l + c^iArlOO-4(— 15)] , (13) 

where (Jq = 2.6 /las, (T15 = 10.2 /xas, and crn = 0.012 describe the systematic-limited, photon-noisc- 
limited, and readout-noise-limited regimes for GAIA, respectively. Gould & Salim (2002) use R 
as a representative magnitude. We apply the same formula to our data in V, and assume their 
parameters without any modification. 

For N observations with uniform sampling conducted over a duration oi tE, the proper motion 
accuracy can be expressed in terms of cxpar as 

= 2V3Ar-VVpar/tE. (14) 

For the GAIA mission, it is predicted that on average N = 82 and Ie = 5 years (Zwitter & Munari 
2003). The expected proper motion errors obtained by combining equations (13) and (14) are 
presented in Figure 3. The sampling will not be uniform for GAIA mission but formula (14) is a 
reasonable zeroth order approximation. 

4. Data 

There have been a few quasar-galaxy associations discussed individually in literature (Burbidge 
1997, 1999; Burbidge & Burbidge 1997; Chu et al. 1998; Arp 1999). We have chosen a different 
approach and make a derivation of proper motions of a large number (8382) of close quasar- 
galaxy pairs from the catalog by Bukhmastova (2001). This catalog is based on the Lyon-Meudon 
Extragalactic Database^ that contains almost eighty thousand galaxies as well as the 8th edition 
of the Quasars and Active Galactic Nuclei catalog by Veron-Cetty k. Veron (1998) that contains 
11358 objects. The associations were selected requesting: (1) the availability of sky positions and 
redshifts for the quasars and galaxies, (2) the galaxy's redshift zq > 0.0004, (3) redshift of QSO 
greater than the redshift of the galaxy {zq > zq), and (4) the projected distance between the 
galaxy and QSO at the galaxy's redshift < 150 kpc for Hq = 60 km/s/Mpc. The Bukhmastova 
(2001) catalog was our only source of data throughout the analysis. 



®http : //leda . imiv-lyonl . f r/sample . html 
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5. Results 

Using the theoretical results from §2, we have constructed a new catalog which evaluates proper 
motions of quasars in QSO-galaxy associations. The complete version of this catalog, which will be 
published in the electronic edition of the Journal, contains 8382 lines and is composed of 9 columns. 
We present 20 randomly chosen entries from this catalog in the printed version of Table 1. Columns 
1 and 4 give the names of the QSO and galaxy, respectively, and columns 2 and 5 the redshifts of 
the QSO and galaxy, respectively. In column 3 we give the visual magnitude Vq of the QSO. These 
were taken directly from the catalog by Bukhmastova (2001). The next four columns contain the 
quantities obtained by us. In column 6 we list coefficients k [Eq. (8)]. These are followed by the 
maximum proper motions for the local scenario jUiocai,max [Eq. (12)] in column 7. We list /Ucosmo,max 
[Eq. (7)] in the 8th column. Finally, the expected proper motion accuracy for GAIA mission, cr^ 
[Eq. (14)], is presented in the last column. This value depends on the magnitude of a quasar, and 
should be also representative of the achievable SIM error. The expected GAIA (achievable SIM) 
accuracy is given for 7126 entries in Table 1. The remaining 1256 objects lack accuracy estimates 
due to one of the following reasons: either Vq > 20.0 (1050 cases), which is fainter than the formal 
limit of astrometric missions, or Vq was unknown (206 cases). 

In Fig. 4 we present the distribution of the number of quasars with a given maximum proper 
motion. The left panel refers to the cosmological scenario and the right panel to the local ejection 
scenario. Typical proper motion in the cosmological case arc <0.1 /ias/yr and the distribution 
of values is rather narrow. This is a consequence of the fact that the number of quasars peaks at 
z ~ 2. The distribution of expected proper motions in the ejection scenario is wider with the proper 
motions expected to be typically larger than a few milli arcsec/yr and reaching 0.1 arcsec/yr in the 
most extreme cases. 

Figure 5 illustrates the ratio of the proper motions in the cosmological and local scenarios on a 
quasar-by-quasar basis. Typical quasars have proper motions in the standard cosmological scenario 
that are 5-6 orders of magnitude smaller than those in the local scenario. This result confirms our 
simple expectation given in eq. (1). 

In Fig. 6, we plot the number of associations as a function of the ratio of the proper motion and 
astrometric errors. The solid histogram presents the results for the cosmological interpretation of 

redshifts and the dotted histogram for the local ejection scenario. We assume Hq = 70 km/s/Mpc. 
The vertical thick line marks the 5a detection limit. For the local scenario, all 7126 maximum 
proper motions should be detectable at the 5a limit. Even if the proper motions are taken to be 
only one tenth of the maximum value, GAIA should detect 99.7% of them; if they are only one 
hundredth of the maximum value, GAIA should still detect 94.4% of them. Figure 2 shows that for 
a typical «; ~ 3, the proper motion stays within a factor of a few from the maximum proper motion 
for a wide range of angles a. Therefore, even if the objects have randomly distributed ejection angles 
which would shift the histogram toward lower signal-to-noise ratios, the proper motions should be 
detectable in the great majority of instances. In the cosmological case, proper motions are too 
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small to be detected with present or presently envisioned technology. The proper motions for only 

10 quasars from our sample (0.14%) have a chance to be detected at their cosmological distances. 
Nevertheless, the non-dctcctions of proper motions would undermine the kinematic interpretation 
of the redshifts of quasars in QSO-galaxy associations, giving support to the complete universality 
of the cosmological interpretation. 

We warn that in some cases the apparent proper motion of quasars may be associated with 
their variability rather than relativistic motion of an ejected object. Such "false positives" could 
imply the local scenario even if the cosmological scenario is correct. Variability of a quasar is the 
sum of the variations coming from the continuum-producing region, broad-line region, and narrow- 
line region. The typical sizes of these regions are 10~^ pc, 0.01 pc, and 10 pc, respectively (Peterson 
1997). If we assume that the sources of the quasar's light can move by the characteristic sizes of the 
above mentioned regions in a coherent way, then the centroid of light from a given region can move 
by about 10~^, 1, and 10^ /xarcsec, respectively for a quasar at a distance of 1 Gpc. Therefore, 
any substantial contribution to flux variability from the narrow-line region may produce spurious 
detection of the quasar proper motion. 



6. Conclusion 

We have designed a proper motion test to distinguish between the cosmological and Dopplcr 
(local ejection) origin of quasar redshifts. We considered the process of relativistic ejection and 
derived a new expression for the maximum proper motion given the rcdshift of a galaxy zq and 
the redshift of an ejected object (quasar) zq. This maximum proper motion depends only on the 
redshifts zq and zq and the assumed cosmological model. 

Using the Bukhmastova (2001) catalog of 8382 QSO-galaxy pair associations we estimated 
typical proper motions for the two considered scenarios. In the standard interpretation of quasar 
redshifts, their typical proper motions are a fraction of micro arc-second, and beyond the reach of 
planned astrometric missions like GAIA and SIM. On the other hand, the quasars ejected from local 
AGNs at velocities close to the speed of light would have proper motions 5-6 orders of magnitude 
larger, and so easily measurable with future astrometric missions, or even in some cases with 
HST, VLT and Keck telescope^. The distributions of proper motions for the cosmological and 
local scenarios are very well separated. Moreover, the division corresponds nicely to the expected 
accuracy from GAIA and SIM. 

We thank Greg Rudnick for a very careful reading of the original version of this manuscript 



The measurements with HST, VLT and Keck telescope would be probably possible only for the brightest 
quasars. The HST Fine Guidance Sensor reaches a per-observation precision of 1 milli arcsec for V < 16.8 
(http://www.stsci.edu/instruments/fgs), and NAOS/CONICA at VLT or NIRC at the Keck telescope reach 
a few milli arcsec precision from the diffraction limited photometry (e.g., Genzel et al. 2003; Ghez et al. 2003). 
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Fig. 1. — The fraction Dm/Dh as a function of redshift for A cosmology with = 0.3 and 
Qa = 0.7. 
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Fig. 2. — The a-dependence of F{a, n) for /t = 2, 3, 4. 
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Fig. 3. — GAIA proper motion errors under the assumption of uniform sampling, the total number 
of observations equal to 82 and the duration of the mission of 5 years. 
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Fig. 4. — The distribution of the maximum proper motions for the quasars from the Bukhmastova 
(2001) catalog in two scenarios for the origin of quasar redshifts. 




Fig. 5. — The distribution of the ratio of quasar proper motions expected from two scenarios for 
the origin of quasar redshifts. 
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Fig. 6. — The ratios of proper motions to the errors of proper motion measurements expected 
from the GAIA astrometric mission. For this plot, we assumed h = 0.7. The sohd histogram is 
based on the cosmological scenario and the dotted one on the local ejection scenario. The GAIA 
5(T-detectability limit is marked with a vertical line. 
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Table 1. Proper Motions for Quasars in QSO-Galaxy Associations 



Quascir 






Galaxy 


ZG 


K 


/^local,max 

[10~®/i~-'arcsec/yr] 


Mco8mo,max 

[ 1 ^ ^ /i ~ ^ arcsec/yr] 


f^j.GAIA 
[10~^arcsec/yr] 


MS 23574-3520 


0.5080 


17.60 


PGC0000621 


0.0007 


1.507 


22707 


0.157 


2.76 


PKS 2357-326 


1.2750 


18.70 


PGC0073049 


0.0008 


2.273 


29971 


0.077 


5.01 


TEX 2358-1-189 


3.1000 


20.50 


LEDA0138064 


0.0024 


4.090 


18758 


0.047 


too faint 


Q 2359-397 


2.0300 


19.00 


PGC0001014 


0.0004 


3.029 


79867 


0.058 


5.99 


Q 0000-398 


2.8270 


18.80 


PGC0001014 


0.0004 


3.825 


100876 


0.049 


5.31 


Q 0000-4244 


1.7000 


20.40 


PGC0001014 


0.0004 


2.699 


71169 


0.064 


too faint 


Q 0000-4239 


2.1900 


21.10 


PGC0001014 


0.0004 


3.189 


84085 


0.055 


too faint 


Q 0001-4227 


2.2400 


19.20 


PGC0001014 


0.0004 


3.239 


85403 


0.055 


6.78 


Q 0001-4256 


2.0300 


18.80 


PGC0001014 


0.0004 


3.029 


79867 


0.058 


5.31 


Q 0001-4225 


1.3000 


20.50 


PGC0001014 


0.0004 


2.299 


60625 


0.076 


too faint 


Q 0001-4255 


2.0400 


18.40 


PGC0001014 


0.0004 


3.039 


80131 


0.058 


4.22 


Q 0002-387 


2.2300 


19.90 


PGC0001014 


0.0004 


3.229 


85139 


0.055 


10.79 


Q 0002-422 


2.7580 


17.21 


PGC0001014 


0.0004 


3.756 


99057 


0.049 


2.27 


Q 0002-4305 


2.2000 


19.80 


PGC0001014 


0.0004 


3.199 


84349 


0.055 


10.06 


87GB 00250+4458 


0.9710 




PGC0001777 


0.0006 


1.970 


41554 


0.093 




Q 0025-4047 


2.1800 


17.14 


PGC0001014 


0.0004 


3.179 


83821 


0.055 


2.19 


Q 0025-4026 


1.1730 


19.00 


PGC0001014 


0.0004 


2.172 


57278 


0.081 


5.99 



Note. — 



The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 



